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Introduction 



When studying evolutionary equations of mathematical physics, we are frequently con- 
fronted with the issue of heterogeneous materials that result in variable coefficients in the 
(partial differential) equations. If the heterogeneities occur only on a small scale, these 
coefficients might be highly oscillatory resulting in the need for significant computational 
effort to solve these equations. That is why one is interested in the behavior of the so- 
lutions when the (small) scale tends to zero. Moreover, if the solutions converge in some 
appropriate sense, it is of interest to know whether the limit is a solution of a similar 
equation and to determine the coefficients. The main objective in homogenization theory 
is to show the convergence of the solutions as the small scale tends to zero and to derive 
the limit equation. If one assumes periodicity in the coefficients, many results are avail- 
able for particular equations, see e.g. |H [HJ |26| as general references. In the non-periodic 
case one cannot expect a similar behavior as very simple equations show, see e.g. |29[ p. 
x, equation (*)]. In this note, we discuss a general compactness result, proving that the 
existence of a limit equation (at least for subsequences) is independent of any periodic 
behavior of the coefficients. Further, we show that this compactness result may be easily 
applied to coupled systems that arise in the area of so-called multi-physics and to equations 
with memory. Apart from that our formulation of the problem is somewhat different from 
the classical approaches to homogenization theory. Starting out with a formulation of an 
equation of mathematical physics being first order in time and space, we deduce a different 
representation of the limiting equations, see also Section [T] for a detailed discussion. One 
may prefer this strategy if there is no appropriate second order formulation. 

In the literature, there are many techniques available that allow the study of homoge- 
nization in the non-periodic case. We mention here the method of .ff-convergence in the 
sense of IT5] 126] or [8] Definition 13.3], which is well-suited for elliptic equations. The 
method of T-convergence is tailored for variational integrals and optimization problems re- 
lated to them, see e.g. [5l[6l[ll]. G-convergence and two-scale-convergence, see respectively 
[HUMUM] an d [21 EE EH], are more general concepts since they can be applied to many equa- 
tions of mathematical physics. Whereas the notion of G-convergence in the very general 
sense of [36] p. 74] (see also Section 2] in this paper) may lead to results that are not precise 
enough, the notion of two-scale convergence may be too restrictive to give useful results 
for the general situation. In particular, n-scale convergence needs additional consideration, 
cf. [T2] p. 2] or [I]. For a possible way of dealing with specific non-periodic coefficients, we 
refer to the generalization of two-scale convergence in |17 | 118 ] . We give a structural way of 
discussing homogenization problems, which was introduced in [291 EQ] with extensions in 
|31| . The idea is to use the observation of |21[ 122] to formulate homogenization problems 
in an operator-theoretic language. 

The core of [2T] and |2"2"] Chapter 6] is that many linear evolutionary problems in mathe- 
matical physics may be written as 



d w + Au = f, 
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Introduction 



completed by some initial conditions and a constitutive relation or material law Ai, being 
a continuous linear operator linking w and u via the equation 

Mu = w. (0.2) 

Here do is a realization of the time-derivative as a normal continuously invertible operator, 
A is a skew-self adjoint operator in some Hilbert space H modelling the spatial derivatives 
and / is a given right-hand side. We assume here that Ai may be represented as a 
function of Bq 1 . This leads to time-translation- invariant and causal material laws. We 
may substitute the constitutive relation (10. 2p into (10.11) to give 

d Mu + Au = f. (0.3) 

This way of considering partial differential equations may be unusual. However, this point 
of view gives a convenient way to discuss many problems in homogenization theory in a 
unified manner. Indeed, the material coefficients M. consist of the conductivity matrix if 
one dicusses the heat equation or the electric permittivity, the conductivity and the mag- 
netic permeability in case of Maxwell's equations. Thus, the discussion of homogenization 
problems in this setting boils down to the discussion of stability under perturbations of M. 
in a suitable topology in equation ( 10.31) . 

Our model of homogenization theory may now be stated as follows, see also Section [T] for a 
more detailed discussion. Consider a sequence of material laws (M. n ) n and corresponding 
solutions (u n ) n of the equation 

d M n u n + Au n = f. 

We ask, whether the sequence (u n ) n converges to some limit v and whether there is a 
material law A/", such that v solves 

d Afv + Av = f. 

In |3T], this question was satisfactorily answered in the case of A having compact resolvent, 
which is the case if considering the heat equation or the wave equation on a bounded open 
set Q satisfying some standard geometric requirements. In addition, it is assumed that 
the material law has a special structure such that the respective equation may also be 
written as a particular second order system. With this approach it is possible to consider 
fractional time-derivatives or ordinary differential equations as consititutive relations, see 
|314 Theorem 4.3 and Theorem 4.5]. In this paper, we generalize these results in the sense 
that we only require A to have compact resolvent when restricted to a domain orthogonal 
to the nullspace of A, i.e., instead of assuming that0 (D(A), r > r > (H, \-\h), we only 
assume the nullspace- compactness-property (or (NC) -property for short) see Corollary 14.71 

1 For Hilbert spaces Hi,Hi and a linear operator A : D(A) E H\ — * H^ with domain D(A), we denote 
the norm in the Hilbert space Hi by [ - 1 and the graph norm of A by | - 1 ^ . If Hi is continuously 
embedded in H2, we write Hi H2 or (Hi, \ -\hj_) ^ (H2, \ '\h 2 )- If this embedding is compact we 
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below for the definition. A relevant example for an operator satisfying the (iVC)-property 

/ curl\ 

is the Maxwell operator A = ( ^ J in L2(£l) 6 , where curio is the curl operator 

with electric boundary condition and Q £ I 3 is a bounded domain satisfying additional 
geometric requirements, cf. [551 Theorem 2] or [T^l [23] • Furthermore, in [2H] it is shown 

that the operator ^ ^ ^Oy defined ^ n L 2 (fl) N+1 satisfies the (iVC)-property, where 

Q ii M. N is open and bounded, grad is the distributional gradient in L 2 (f2) with domain 
equal to W r 2 1 (fi). The index "0" refers to Dirichlet boundary conditions and div is the 
negative adjoint of grad . The same reasoning can be applied to the spatial operator of the 

elastic equations ( ^ Grad \ ^ w j iere Q ra( j j s ^ e symmetrized gradient as defined in 



v - Grad 

Definition 4.9] with some boundary conditions on a bounded domain Q satisfying suit- 
able geometric requirements, cf. e.g. [33"1 Theorem 2]. Thus, our main theorem, Theorem 
14.51 ma Y be seen as a general theorem giving a compactness result for the homogenization 
of (coupled) equations in mathematical physics. We shall also mention that the results 
obtained in this article not only generalize [29| Theorem 2.3.14] but improve the represen- 
tation of the homogenized equations. Moreover, we show that the homogenized equations 
satisfy the assumptions of Theorem 12.11 i- e -> we have a solution theory for these equations 
and causality of the solution operator. A detailed discussion of our main result is given in 
Section [TJ In this section we also give an account of the ideas used and compare it to other 
strategies in the literature. 

We sketch our plan to achieve our main result as follows. In Section [2J we discuss the 
mathematical framework of evolutionary equations and recall the main theorem of |21] . 
Section [3] sketches the ideas, definitions and main theorems of RU] and 1311 . In Section [4] 



we present our main result, Theorem 14.51 We show optimality of this theorem by means 
of counterexamples. The proofs in Section [4] also refer to results from Section [6], where 
some technical lemmas are provided. The results from Section [6] are needed to prove 
well-posedness of the limiting equation constructed in Theorem 14.51 Before, however, 
completing the proof of Theorem 14.51 in Sectional we apply Theorem 14. 51 to some examples 
from mathematical physics in order to illustrate our findings in Section [5l 

We indicate weak convergence in a Hilbert space by '— »' or 'w-lim'. Norm-convergence 
will be denoted by '— >' if not specified differently. The underlying scalar field of any vector 
space discussed here is C. 



1 Discussion of the main result 

Starting out with a different formulation of the problem similar to f l0.3p . we derive ho- 
mogenization results for a large class of partial differential equations. The homogenized 



write H x — + H 2 or (H lt \-\ Hl ) (H 2 , \-\ Ha 
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1 Discussion of the main result 



coefficients, however, have a different representation as in classical results in the literature, 
cf. e.g. jUEl]. We illustrate the difference between the classical approach and the approach 
considered here with the heat equation (cf. e.g. [221 p. 350]): 

Let SI g M. N be a bounded domain. Denote by 9 : (0, go) x Q — > R the temperature 
distribution and by q : (0, go) x SI — > M. N the heat flux. The heat equation is a system of 
the two equations 

d 9 + div q = f 
q = —Kgrad9. 

Here 3q9 denotes the time-derivative of 9, f is a given source term and k: Q — > M. NxN is a 
bounded function being the (material dependent, symmetric) conductivity tensor satisfying 
k(x) — c for some c > and all iefi. The first equation in (11. ip is called the /ieai /Ziw; 
balance and the second one Fourier's law. The system is completed by boundary and 
initial conditions. For simplicity, let us assume Dirichlet boundary conditions for 9 and 
zero initial conditions. Similar to the introduction, as a reminder for Dirichlet boundary 
conditions, we shall write grad instead of grad. Of course, 9 and q are the unknowns in 
the system. 

The classical way of discussing the heat equation is to plug Fourier's law into the heat flux 
balance in order to arrive at 

d 9 -dwKgrad 9 = f. (1.2) 

Now assume that k has an extension to M. N such that k{x + ej) = k(x) for all x e M. , Cj := 
(<%)ie{i,...,Ar}> 3 E {I? • • • ) N}. In other words, assume that k is (0, 1) ^-periodic. In homog- 
enization theory one is now interested in the behavior of solutions coming from equations 
having highly oscillating coefficients. A possible way to model that is to discuss a sequence 
of problems, i.e., for neN consider the solutions (9 n , q n ) and 9 n of the following respective 
equations 

d 9 n + div q n = f 
q n = -K(n-)gmd 9 n . 

and 

d Q 9 n - div n{n-) grad 9 n = f. (1.4) 

Standard a priori estimates imply that (possibly after passing to a subsequence) (9 n , q n ) n 
and (9 n ) n converge weakly to some functions (9, q) and 9, respectively. In order to de- 
termine 9, classical results derive the temperature distribution to be the solution of the 
equation 8q9 — div«;og ra do# = / with the homogenized (constant-coefficient-) matrix k 
from equation (II. 4p . A main difficulty to overcome in that approach is the question of how 
to perform the limit in J n k(jix) grad 9 n (t, x) ■ grad 9 n (t, x) dx. A possible way to deduce 
convergence of the integral expression to $ Q k,q grad 9(t, x) ■ grad 9{t, x) dx is the famous 
"div-curl-lemma" or the "compensated compactness", which is due to Murat and Tartar, 
see e.g. [T3J 125]; or |2SJ U\- The strategy of doing so is well-established and can also be 
applied to other cases such as linearized elasticity. 
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Starting out with the sequence of equations given by (jl.3p . we propose another way of 
deducing the limiting equation. In the following lines we only sketch the ideas as they are 
rigorously proven in the forthcoming sections. Written in a block-operator-matrix-form 
the equations (11. 3ft read as 

<;K<M + u?))ft)-(o- <-> 

Following the introduction and recalling that 8q can be realized as a continuously invert- 
ible operator, the latter equations are clearly of the form (10. 3 p with Ai n = + 

3q 1 \ n ( J and A = [ ^ ^ A ^J. Note that, due to the Dirichlet boundary con- 



v «(n-) J ygradg 

ditions, A is a skew-selfadjoint operator in L 2 (Q) N+1 since div* = — grad . Having a 
homogenization result for A having compact resolvent at hand (see |29l Theorem 3.5] or 
Theorem 14.31 in this article), we want to deduce our homogenization result from the com- 
pact resolvent case as treated in Theorem 14.31 Due to the infinite-dimensional nullspace 

/ div\ 

of div in spatial dimensions greater than 1, the operator , never has a com- 



,grad 

pact resolvent if N ^ 2. But the domain of grad equals W 20 (Q), which is, if endowed 
with the graph norm, compactly embedded into L 2 (Q) due to the selection theorem by 
Rellich and Kondrachov (keep in mind that fl is assumed to be bounded). In that sense 
the operator grad is the 'good' part, whereas div is the 'bad' one. To overcome this 
problem the idea is to restrict A to a domain being orthogonal to the nullspace of A. 
Due to Dirichlet boundary conditions the operator grad is one-to-one, thus N(A), the 
nullspace of A, equals {0} iV(div) E L 2 {Q) ® L 2 (Q) N . Since div* = — grad , we have 
L 2 (Q) N = P(grad ) © iV(div). Using Poincare's inequality, we deduce that the range of 
grad is closed in L 2 (Q) N . Hence, 

L 2 (Q) N = P(grad )©iV(div). 

Along that decomposition of L2(£l) N , we decompose the heat flux q n from equation (II. 5p as 
Qn = Qn + Qn with q\ e i?(grad ) and q\ e iV(div). A convenient way of writing the operator 
A acting on the three components (0 n , q n , q%) is by introducing the operator P: L 2 (fl) N — > 
i?(grad ), which maps any g e L 2 (fl) N to its orthogonal projection Pg in the range of 
grad . Looking at P in that way, the adjoint of P is the identity with domain i?(grad ) 
and target space L 2 (Q) N . Hence, A as an operator acting on L 2 (Q) © i?(grad ) © iV(div) 
may be written as follows 

/ div P* N 
A = Pgrad 
\ 0, 

Observe that A leaves the space L 2 (Q) © P(grad ) invariant. Moreover, note that A is 
one-to-one and skew-selfadjoint on L 2 (Q) © P(grad ). Furthermore, it is not hard too see 
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1 Discussion of the main result 



that the domain of divP* if endowed with the graph norm is compactly embedded into 
P(grad ), see e.g. |3T] Lemma 4.1]. Thus, the operator L ^ ^ J has compact 



resolvent. 



Denoting by Q: L>2(fl) N — >■ iV(div) the operator, which maps g e L2(fl) N to its orthogonal 
projection Qg e iV(div), we deduce from (jl.5p the following system 

oo o \ / o divP*o\\[e n \ ft 

Pn{n-)- x P* Pft(n-) _1 Q* I + I Pgrad ] | ( q l n \ = ( 





^0 0^ 




N 









{0 0) 





K Q^n^P* QK(n-)- l Q*J \ / / \&J \0, 

(1.6) 

Now, we could apply Theorem 14.31 to the first two rows of equation (11. 6p and write the 
term PK(n-)~ l Q*q 2 l to the right-hand side, i.e., 

l o\ /o o W o div p*\ \ / e n \ ( f 



0) + lo PK(n-)- 1 P*J + LPgrad J J V -PK{n-)- x Q*q 2 



A suitable choice of subsequences would then yield the weak convergence of Pn(n-)~ l Q*q 2 l . 
In order to get a similar equation as in the original system after letting n — > oo, we need 
the limit of P n(n-)~ l Q* q 2 , to be of the form Cq 2 for a suitable bounded linear operator C 
and q 2 being the limit of {q 2 )n- The reason for that is that the limiting equation should 
contain the limit of q\ as an unknown. However, in general, we cannot deduce that the 
limit of P ^{n-)^ 1 Q* q^ is the product of the limits of PK,(n-)~ 1 Q* and q 2 . A deeper look 

divP 



into Theorem 14. 31 reveals that due to the compactness of the resolvent of . 

1 Pgrad 

the sequence (8 n , ql) n converges in a stronger sense than only weakly such that a particular 
'weak-strong principle' can be applied, cf. Theorem 13.51 Hence, we want to express q 2 in 
terms of (9 n ,q^). Therefore we are led to perform Gauss steps to transform ( II. 6p into a 
system such that the convergence of (q 2 ) n can be deduced from the convergence of (9 n , q^) n - 
A possible way of doing so is the following 



~ X Q* (QKin-)-^*)- 1 Q K (n- 
Qk(u-)- 1 P* 









(0 







H 


PK(n-) _1 P* - Pk{u 




\0 0) 




VO 



divP* N 
| Pgrad I I I I = I 
0// \q 2 J Vo, 




'10 
Multiplication by | 1 ] gives 

,0 (QK(n-) _1 Q*)" 
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^10 0^ 











h 


P«(n-) _1 .P* 




Vo o o) 




Vo 




N 
- PK(n-) _1 g* (Q/^n-^QT 1 Q«(n-) _1 P* 

o di V p* o\\ fe n \ (t 

In the latter equation we may perform the limit at least for a suitable choice of subse- 
quences, for which we use the same notation. To this end, note that the operators 

'0 0\ 

Pk(u-)- 1 P* - P^n-^Q* (Q^n-) -1 ^*) - * QK{n-)- l P* (n £ N) 
,0 (QK(n-)- 1 Q*)~ 1 QK(n-)- 1 P* 1/ 

form a bounded sequence in the space of linear operators in the separable Hilbert space 
L2(Q) N+1 . Thus, there exists a subsequence, which converges in the weak operator topology 
of L (L 2 (fi) Af+1 ). Applying Theorem 14.31 to the first two rows of the latter equation, i.e., 



i o Wo o 

OJ \p P^n-^P* - Puin^Q* (QK(ri-) _1 Q*) Qn{n-)- l P 

divP*\\ (6 n \ 



+ 



Pgrad J J \q}J V 



we deduce that (6 n , q^) n converges weakly and due to the compact resolvent of the oper- 
divP*\ + , ^ f ._ 3Q ._ 3 , 



ator ^pg ra( ^ g J ^e sequence (d 9 n (t),d converges strongly in L 2 (£l) © 

i?(grad ) for all t £ M, cf. Theorem 14.31 In particular, this means that 

{[QK{n-)- 1 Q*Y 1 QK{n-)- 1 P*q]) n 

converges to the product of the limits of (^(Q K,(n-)~ l Q*) 1 Q K,(n-)~ l P*\ and (<7*) n > see 
also Corollary 13.61 below. Hence, in the last row equation 

{QKin-y'QT 1 Q<n-)- l P*ql + & = 

it is possible to let n tend to infinity. We get that (0 n , q^, q£) n weakly converges to a 
solution of the following equation 



1 N 
d | 
0, 



/0 0\ 

lim^oo (Pk(w)- 1 ?* -P«;(nO _1 Q* (QK(n-) _1 Q*) _1 QK(n-) _1 P*) 
v \^ n ^ [{Q Kin-)- 1 Q*Y 1 QK{n-)- l P*^ 1 J 



11 



1 Discussion of the main result 




Note that the system is of a rather simple form, since the right-hand side contains two 
zeros. In general, we have a more difficult situation, which then also results in a more 
involved limiting equation. 

How is it possible to see that the limiting equation is well-posed? This requires the per- 
formance of other Gauss steps and so the well-posedness is not easy to see. Following the 
strategy of Section [Gj, we end up with the following equations 



- f 1 o o\ 
d o o o 
\o o o/ 



• 

+ 



lim n ^ao(PK(n-) _1 ^*--PK(™-) _1 Q*(QK(n-)- 1 Q*) ^(n-)"^*) limbec (P^n-^Q* (Q/^n-)"^*) 
lim n ^. X: ((QK(n-)- 1 Q*) _1 QK(n-)- 1 P*) lim^a, (QK(n-) _1 Q* ) ~* 



+ ('^T;)) (?)-«)■ ^ 

The latter equation can be shown to be well-posed, time-translation invariant and to have 
a causal solution operator. Note that, if one is only interested in the behavior of the 
temperature distribution, we can reformulate the latter equation into a second order form. 
The resulting equation would be 

d 9-divP* ( lim PK(n-) -1 P* - Pn(n-)- x Q* (QK(n-)- 1 Q*Y 1 Q«(n-) _1 P* )~ Pgrad 9 = f. 

Vn— >oo v ' J 

Using the periodicity of k, we deduce that Pn(n-)~ l P* converges in the weak operator 
topology to Ph Q1 \N ^(x)^ 1 dxP*, cf. [3T| Proposition 4.3]. Thus, in order to answer the 
question proposed in the title of this article, the (inverse of the) classical homogenized 
matrix and the harmonic mean of k differ on the subspace i?(grad ) by the limit of the 

sequence (^Pn(n-)~ 1 Q* {Qk^u^^Q*) 1 Qn(n-)~ 1 P*^j in the weak operator topology. 

The interested reader might think, why such a seemingly complicated strategy yielding the 
homogenized equations should be applied. In the case of the heat equation this strategy 
indeed does not give anything new despite the fact that the homogenized equations have 
a different representation. For the heat equation it is even worse: with this strategy one 
cannot deduce the convergence of the whole sequence. However, note that the approach 
presented here only uses abstract theory from functional analysis and does not rely on the 
specific form of k being a periodic multiplication operator. If k is a linear operator invoking 
non-local terms, well-known homogenization theory might fail to work. Moreover, the 
way of computing the homogenized coefficients carries over to a large class of evolutionary 
equations: It is possible to treat Maxwell's equations, the wave equation, the heat equation 
or general coupled systems in mathematical physics in a unified manner. It is also possible 
that a second order formulation might not be available or is not easy to handle, cf. e.g. 
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[2"2"| Equation (6.3.9), p. 455] so that the usual strategy might not work. 

Now, we give a detailed account of the rigorous treatment of evolutionary equations in the 
abstract setting of equations ( 10. ip and ( 10. 2p . 

2 Setting 

We recall the setting of evolutionary equations established in [21] or |2"2"j Chapter 6]. Let 
H be a Hilbert space and denote by L 2 (R; H) the space of H- valued L 2 -functions. The 
operator 

d : W^(R; H) g L 2 (R; H) L 2 (R; H) : f ^ f 

assigning to each weakly differentiable iJ-valued function its weak derivative is skew- 
selfadjoint. Define the unitary Fourier transform J 7 : L 2 (R: H) — > L 2 (R; H) as the closure 
of the mapping 




defined for / belonging to C™(R;H) the set of indefinitely differentiable, if- valued func- 
tions with compact support. Let 

m : {/ 6 L 2 (R; H); (x - xf(x)) e L 2 (R; H)} g L 2 (R; H) - L 2 (R; H) : f ~ (x ~ 

For i/ > define ff^o(M; if) := L 2 (M, exp(— 2vx) dx; H) the space of //-valued (equivalence 
classes of) square-integrable functions with respect to the weighted Lebesgue measure with 
Radon-Nikodym derivative exp(— 2z/(-)). We also write H V $(R) if H = C. The mapping 
exp(— vm) : H u q(R; H) — > L 2 (R; H) : f (x h- » exp(— ux)f(x)) is unitary and the 
operator 

d 0i j, : = exp(— vm)*(d + v) exp(— vm) 

is normal in H U ^(R;H). If there is no risk of confusion, we simply write do instead of 
do,v We have 3qI e L(H Ut o(R; H)) with = 1/v. Introducing the Fourier- Laplace 

transform C u : = J 7 exp(-vm), we get 

do,v = C* v (im + v)C v . 

Consequently, 

The latter equation gives a functional calculus for the normal operator d^l'. 

Definition (Hardy space and functional calculus for do tU ). For an open set E g C and a 
Banach space X, we define the Hardy space 

H <X '{E] X) := {M : E -> X; M bounded, analytic} 
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2 Setting 



and \\M\\ X := sup{\M(z)\ x ; z e E). Let H u H 2 be Hilbert spaces, v > 0, r > l/(2i/). 
• For M e W°(B c (r, r); L(H U H 2 )) define 



M < a »>= £ *( M fe)) £ "' 



where (M 0) (() » M (g^) 0(i) for all t e R and <j, 6 Cf (K; #,). 

For c > define 

TT'XBc^r);^)) := {M e n m (B c {r, r); L(^)); Rez" 1 M(z) ^ c (2 e B c (r, r))}. 



For easy reference, we call elements of 'H oc (i?c( r ) r )', L(Hi)) material laws or constitutive 
relations and elements of ^°°' c (i?c(^, ?"); L(H\)) (c) -material laws. 



As the definition of do indicates, we will deal with space-time operators. In particular, we 
identify any closed, densely defined operator A : D(A) E H — > H in some Hilbert space if 
with its canonical extension on the space of if- valued H u (lR)-functions, cf. |2Tj. We have 
the following well-posedness theorem taken from [21J. 

Theorem 2.1 ([211 Solution theory]). Let H be a Hilbert space, c, v > 0, r > l/(2v) and 
M g W°> c (B(r, r); L(H)). Let A : -D(A) ^ H ^ H be skew-selfadjoint. Then the equation 

{d M{d l ) + A)u = f 

admits a unique solution u e H Ut0 (W; H) for a dense set of right-hand sides f e H Uj q(M.; H). 
Moreover, the solution operator (5oM(5 _1 ) + A)~ 1 is a densely defined, continuous operator 
in H U fi(R; H) with operator norm bounded by -, and the operator (5oM(5 _1 ) + A)^ 1 is 
causal, i.e., for all f e H U fi(R; H) and a e R we have 

X(-o O ,a)(m)(a M(4 1 ) + A)-\ X (-^a){m)f) = X ^a){m){d,M{d Q l ) + A)^(f), 

where X(-oo,a)("^) denotes the multiplication operator mapping f e H Ut0 (M; H) to the trun- 
cated function 1 1— > X(-x<,a)(t)f(t). 

We note that these results carry over to "tailor made" distribution spaces - so-called Sobolev 
lattices - discussed in (20]. In |31[ Remarks 1.2: (i)-(iii)] and |21[ Sections 2 and 3] the 
core issues are sketched. We will use the notation from [21] and for the sake of clarity, we 
recall the main definitions. For k 6 Z, a Hilbert space H and a densely defined, closed 
linear operator C : D(C) E H — * H with e g(C), we denote by H k (C) the Hilbert space 
defined as the completion of D(C' fc ') with respect to the norm |-|# fc (c) : u l— * \C k u\u- It can 
be shown that the closure of H\ k \(C) E H k {C) — * H k _i(C) : u >-> Cu is unitary. We will 
re-use the name C for this extension. We are interested in the special cases C = A + 1 with 
A skew-selfadjoint or C = 8q. For £ e { — 1,0, 1} we let H^a := H^A + 1). For 8q defined 
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on i/ I/i o(M)-functions with values in a Hilbert space H we write H U ^(R; H) := H k (do). 
Consequently, we also use the spaces H^k(R; H^a), f- e { — 1, 0, 1}. The extension of the 
solution operator to H v _i(R;H) also serves as a way to model initial value problems, see 
e.g. |22| 6.2.5 Initial Value Problems]. 

3 Preliminary results 

We summarize some findings from [3T)| I3"T] . 

Definition. For an open set E E C and Hilbert spaces Hi,H 2 , we define on the set 
r hL co {E; L(H\, H 2 )) the initial topology t» induced by the mappings 

U*{E- L(H U H 2 )) b M ~ <0, M(-)V> e n(E), 

where H(E) is the set of holomorphic C-valued functions endowed with the compact open 
topology. We define H™(E; L(H U H 2 )) := {W°{E;L{Hi,H 2 )),tm) and re-use the name 
W%(E; L(Hi, H 2 )) for the underlying set. 

Theorem 3.1 (sequential compactness, [30| Theorem 3.4]). Let Hi, H 2 be separable Hilbert 
spaces, SiC open. Let B g H™(E; L(Hi, H 2 )) be bounded, i.e., sup{||M(z)|| ;z<=E,M<e 
B} < oo. Then B is relatively sequentially compact. 

Lemma 3.2 ([201 Lemma 3.5]). Let H be a Hilbert space, r > 0. Let (M n ) n be a 
bounded and convergent sequence in the space H^(B(r, r); L(Hi, H 2 )) with limit M e 
H™(B(r, r); L(Hi, H 2 )). Then (M re (<5^~ )) n converges to M(3q 1 ) in the weak operator topo- 
logy of L(H Uik (R; Hi), H Vfk (R; H 2 )), where v > l/(2r), fceZ. 

Proof. In [301 Lemma 3.5], the claim was shown for the case k = and Hi = H 2 . The 
general case follows by observing that <9q : H Ut k(R;Hi) — > H u _o(R; Hi) is unitary and 
obvious modifications. □ 

Lemma 3.3 ( |29[ Lemma 1.5]). Let Hi,H 2 be Hilbert spaces. Let E E C be an open 
disc with centre z and let (M n ) n = (Yuk=o(' ~ z ) k A n k)n be a convergent sequence in 
W™(E; L(Hi, H 2 )) with limit 5]*°=o(' — z ) kj ^k- Then A nk — > A k as n — > oo in the weak 
operator topology of L(Hi, H 2 ) for all k e N . 

For a Hilbert space H and v > 0, we define 

C U (R; H) := {(f) e C(R; H); swp\exp(-vt)(/)(t)\ H < oo}. 

We endow C U (K;H) with the norm \'\c v '■ <P l— * su PteRl ex P( — vt)(f)(t)\H. Recall from |2"2"| 
Lemma 3.1.59] that H Ut i(R; H) is continuously embedded in C U (R; H). 
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3 Preliminary results 

Lemma 3.4 ([211 Lemma 2.2]). Let H be a Hilbert space, v > 0. // (f n ) n in H Uj i(M.; H) is 
bounded and converges pointwise to some f e H Ut i(ti&] H), then 

do'Ut) 5 -V(*), 

/or a// t e K. 

Theorem 3.5 (weak-strong principle, [3T| Theorem 2.3]). Let H be a Hilbert space, e > 
0, (M n ) n be a convergent sequence in 7{^(Bc(0, e); L(H\, H 2 )) with limit M. Then, for 
v > 2/e and any bounded sequence (v n ) n in H u i(R; Hi) and v e H Vt x(R; Hx) such that 
v n (t) v(t) in Hi for all t e W, 

w- \im(M n (d^)v n )(t) = (M(<V»(t) e H 2 , 

n— >oo 

for all t g R. 



Proof. In [31] the proof is given for the case Hi = H 2 . The assertion follows analogously 
with obvious modifications. □ 

Corollary 3.6. Let Hx,H 2 be Hilbert spaces, e > 0, (M n ) n be a convergent sequence in 
U%(B c (0,e);L(H u H 2 )) with limit M e H™(B c (0, e); L(H U H 2 )). Let v > 2/e, k e Z and 
let {y n ) n be bounded in H v ^{R\Hi), v g H u ^(R',Hi). Assume there is I e No such that 
do l v n g H Wi i(R] Hi) and d^ l v n {t) d l v(t) in Hi for all t e R. Then 

w- lim M n {3^)v n = M(d Q l )v g H v>k (R; H 2 ). 



Proof. Since (M n (d 1 )v n ) n is bounded in H u j,{M; H 2 ), there is a subsequence with in- 
dices (nj)j weakly converging to some w g H u ^(R',H 2 ). The assumption guarantees that 
(8q ~ l v n ) n is bounded in H V) i(M.', Hx). Moreover, by Lemma l3.4[ (d^ l v n ) n converges 
pointwise to <9o~' fc ' V Thus, by Theorem 13.51 and the weak continuity of point-evaluation, 
we deduce that, for t el, 

(d lkH w)(t) = w- lim(5 - |fc| -'M^(5 - 1 )^.)(t) 

3— »oo 

= w-\imM nj {d 1 )d~^~ l v n3 {t) 

j->cc 

= M(d^)d W Mt) = d ^ l M{d^)v{t). 
Hence, w = M(5 " 1 )?;. □ 
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4 A general compactness theorem for the 
homogenization of evolutionary equations 



We introduce the concept of G-convergence to bridge the gap between the classical approach 
in homogenization theory and the Hilbert space perspective discussed here. 

Definition (G- convergence, |3SJ p. 74]). Let if be a Hilbert space. Let (A n : D(A n ) g 
H —*■ H) n be a sequence of one-to-one mappings onto H and let B : D(B) g H — > H be 
one-to-one. We say that (A n ) n G-converges to B if for all / 6 H the sequence (A~ 1 (f)) n 
converges weakly to some u, which satisfies u e D(B) and B(u) = f. B is called a G-limit 
of (A n ) n . We say that (A n ) n strongly G-converges to B in H, if for all weakly converging 
sequences (f n ) n in H, (A~ l (f n )) n weakly converges to some u, which satisfies u e D(B) 
and B(u) = w- lining f n . 

Proposition 4.1. The G-limit is uniquely determined. 

Proof. Let H be a Hilbert space. Let (A n ) n be a sequence of one-to-one onto mappings 
which is G-convergent to the one-to-one mapping B : D(B) E H — > H. Define C : = 
{(«, f) e H (& H;u = w- lim^oo A~ 1 (f)}. Then (7 I 5, so that C is a mapping. Moreover, 
since C is onto and B is one-to-one, we conclude that C = B. □ 

Remark 4.2. If (-*4 n ) n in the above definition is in addition a sequence of linear and closed 
operators and B is also closed and linear, then the above definition of G-convergence is 
precisely convergence of the resolvents in the weak operator topology, which is the original 
definition in [36] in the Hilbert space setting. 

We now prove compactness results concerning G-convergence for operators that may be 
associated with evolutionary equations. More precisely, we will deal with the following 

cases: 

Definition. Let Hi,Bi2 be Hilbert spaces. We say a pair ((M n ) n ,A) satisfies 

(PI) if there exists e, r, c > such that (M n ) n is a bounded sequence in 7/ 00 (5(0, e); L(Hi))n 



-H°°' c (5(r,r);L(ifi)) and A : D(A) g ffi — H x is skew-selfadjoint and the embed- 
ding (D(A), \-\a) ^ (Hi, \-\hx) is compact, 



(P2) if there exists e,c,r > such that (M n ) n = ( ( 




is bounded in 



H™(B(Q,e)-L(Hi® H 2 )) nW°> c (B(r,r)-L(Hi® H 2 )) and A = ^ 
((Mii >n ) n , A) satisfies (PI). Moreover, 
(i) for all n e N, R(Mi(0)) = R(M n (0)) and M n (0) ^ c on R(Mi(0)), 



is such that 
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4 A general compactness theorem for the homogenization of evolutionary equations 

(ii) denoting by qj : Hj —*■ R(tt*) n N(M 1 (0)) (j e {1,2}) the canonical ortho- 
projections, we have for all n e N 

(MWO)^)" 1 ^M^ n (0)qtY = qi M[ 2in (0)q* (g 2 M^ n (0)g 2 *)^ . 
With these definitions, the core result in [31] now reads as follows. 

Theorem 4.3 ( |31[ Theorem 3.5]). Let H be a Hilbert space and assume that ((M n ) n , A) 
satisfies (PI) and that (M n ) n converges to N s H™(B(0, e); L(H)). Then there exists 
Vq ^ such that for all v > u , (doM^^Q 1 ) + A)k strongly G-converges to doN^d^ 1 ) + A 
in H v -i(R] H). Moreover, N e H^ c {B{r, r); L(H)) and 

d 3 (d M n (d l ) + A)-\f n (t) - dv 3 (d N(d^) + A)-\w- lim f n )(t) e H 

n— »oo 

ds ti — > oo /or all t e W and all weakly convergent sequences (/„,)„ m i/ I/j _i(R; 5). 

The generalization of this theorem to the case (-P2) requires a homogenization result for the 
case of A = 0, i.e. a result on the homogenization of ordinary integro-differential equations. 
Since we deal with a possibly degenerated case in the sense of [21] 3.3 Some special cases], 
we cannot use the homogenization result for ordinary integro-differential equations already 
established in [30] Theorem 5.2]. The refined argument is tailored for the O-analytic case 
(cf. Section [6]), which, however, does not cover the results in |30j. 

Theorem 4.4. Let H be a separable Hilbert space, e, c,d,r > 0. Let (M n ) n be a bounded 
sequence in "H^i^O, e); L(H)) n T-L aD ' c (B(r,r); L(H)) and assume that for all n e N, 
M n (0) ^ d on i?(M n (0)) = R(M 1 (0)). Then there exists r' e (0,r] and a strictly monotone 
sequence of positive integers (n*.)*, such that, for v > l/(2r'), (c^M^c^" 1 )^ G-converges 
to dofJ-idQ 1 ) in i? l/; _i(R; H), where fi has the following properties: there is e',c' > such 
that 

(i) fi 6 H x (B(0, e'); L(H)) n H^ c '(B(r', r'); 
(«J J R(/i(0)) = i?(M 1 (0)) ; 

(raj /or open I?iC relatively compact in B(0, e')\{0} (: <=^ E cc 5(0, e')\{0}j, 

M^-)" 1 - MO -1 e 7C(£; L(5)) (A - oo). 

Proo/. Define the Hilbert spaces H x := fl(Mi(0)) and H 2 := JV(Mi(0)) together with the 
canonical (orthogonal) projections 7Tj : 5 — > 5j, j e {1,2}. Then, for all n £ N and 
j, 6 {1,2}, set Mjk, n {-) '■= 7r i^n( , ) 7r fc- Now, the first assertion in Lemma 16.101 ensures 
the existence of e' > such that, for all E ii C relatively compact in 5(0, e')\{0}, the 
sequence (M ri (-)~ 1 )„ is bounded in K^iE; L(Hi © H 2 )). By a-compactness of 5(0, 
and Theorem l3.lt we may choose a subsequence (M n)t (-) _1 )fe of (M n (-) _1 ) n such that there 
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is a holomorphic mapping 77 : B{0, e')\{0} — ► L(H) with 

M^-Y 1 - ^ e 7C (^1 )) (* -^»,£cc 5(0, s')\{0}). 

The residue theorem ensures that the residues of M nfe (0 _1 converge to the one of 77 in the 
weak operator topology. Moreover, it is easy to see that the coefficients of the Laurent 
series expansions of M nfc (0 _1 converge in the weak operator topology r w . With the help 
of the first assertion of Lemma I6.10[ the Laurent series expansion of rj is of the form 



(r w -) lim^^ Mn jUk (0) - 1 + M u (z) M 12 (z) _ 

M 21 (z) z-\t w -) lim^oo M 22tnh (0)~ 1 + M 22 (z) 



for suitable bounded holomorphic operator- valued functions Mj^ for j, k e {1,2}. The 
second assertion of Lemma 16.101 yields the existence of e" > such that /i := rj(-) _1 e 
'H co (B(0, e"); L(H)). Moreover, from the representation in Lemma [6.101 we read off that 
R(Mi(0)) = i?(/i(0)) and /i(0) ^ d! on H\ for some d' > according to Inequality 
(16. 3p and the fact that positive definiteness is preserved under limits in the weak op- 
erator topology. Similarly, Re//(0) ^ c' > on id 2 . Thus, by Remark 16.31 it follows 
that yU lies in "H 00 ^ (B(r',r'); L(H)) for some r', c" > 0. It remains to show the G- 
convergence result. To this end let v > l/(2r'). By the convergence of the coefficients 
in the Laurent series of (M nfe (0 _1 )fc, we get that ((0^n*(0 )* converges to (•)?/(•) in 
H™(B(l/(2v),l/(2u))-L(H)). Thus, Lemma K2\ implies that ((doM^d^ 1 )) -1 )*. converges 
to d$ T){dQ ) i n the weak operator topology of L(H„-i(lSL; H)). Employing Remark 14.21 we 
obtain the desired G-convergence. □ 



Theorem 4.5. Let Hi,H 2 be separable Hilbert spaces. Assume that ((M n ) n , A) satisfies 
(P2). Then there exists u > 0,e',c' > and (rik)k a strictly monotone sequence of 
positive integers such that for all v > vq the sequence (<9oM nfe (d^ 1 ) + A)k G-converges 
to (doNido 1 ) + A) in H u -x(R; H x H 2 ) with 



N{.) 



%(0 _1 %(0 V2(-)~ 1 , 

e ft 00 (5(0, £'); L{H X H 2 )) n H x ' c ' (B(l/(2v Q ), l/{2v ))-L{H 1 H 2 )), 

where 

:= lim M n , nfc (.) - M 12jnk (.)M 22jnk (.)- 1 M 2 x > n k (-) e K(B(0, e')- L(#0) 

fe— >oo 

!»(•) := hm (M 2 2,n fc (-)) _1 e U*(E;L(H 2 )) {E S(0,eO\{0}) 

fe— >00 

773(0 := hmM 22infc (-)- 1 M 21infc (-)E^(S(0,e / );^(^i,^2)) and 



fc— >oo 



774(0 := limM^OM^O^e H™{B{0,e'); L{H 2 , Hj). 
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4 A general compactness theorem for the homogenization of evolutionary equations 



Moreover, R(N(0)) = i^M^O)). 



Proof. By Theorem 16.111 (applied to M = M n and the sequence N just the constant 
sequence consisting of M n as every entry) there exist e', r', d > such that, for all neN, 



1 -Mi2,„( 




)M 22 , n ( 
1 



Mu. n 



Mii, n (-) M 12 , n (- 
M 21 ,„(-) M 22 , n (- 

(•) - Mi 2 , n (-)M 22 , n 




1 



-M 22 , n 
_1 M< 



21,nl 



21,nl 



M 22j „( 

e ft°°(5(0, e'); L{H X # 2 )) n W ^ {B{r' , r'); L(fT a # 2 )). 



Let v > l/(2r'). By Theorem 13 . 1 1 and Theorem l4.41 we may choose convergent subsequences 
of the material law sequences 



(/^3,n)n 
(/^4,n)n 



(M n ,„(-) - M 12in (-)^22,„(-) _1 M 21 , n (-)), 

(Maa.nCO -1 )™ 

(M 22in (-)- 1 M 21ire (-)) n and 

(M 12 , n (-)M 22in (-) _1 )n- 



We will use the same index for the subsequences and denote the respective limits by rji, rj 2 , t] 3 
and 774. Using the representation from Theorem 16.61 we get with the help of Theorem 16. Ill 



(Ci e {0}) (G 3 {0}) = R 



( 



\ 



f(M® n 




Mil 





MS?. 



oJJ 



R 
R 



0, 

M u , n (0) M 12 , B (0)^ 
M 2 i,„(0) M 22 , n (0) / 

M n , n (0) - M 12!n (0)M 22>n (0)- 1 M 21) „(0) 

M 22 , n (0) / 

f? (Mu, n (0) - M 12i „(0)Af 22 ,„(0)- 1 M 21in (0)) i2 (M 22 , n (0)) 



Now, Afu,„(0) - Mi 2jn (0)M 22iri (0)^ 1 M 2 i i „(0) is strictly positive on G x = R(M u ,i( )) and 
M22,„(0) is strictly positive on G3 = i?(M 22i i(0)) uniformly in n. Hence, we deduce that 
R(r]i(0)) = H(Mii,i(0)) and, from Theorem Ol that R{r} 2 {-)~ 1 (0)) = fl(M 22j i(0)). Let 
C/1,/2) 6 (R;^ # 2 ) and for ra e N, let ( Uli „, u 2 , n ) e H^R; H x F 2 ) be the 

unique solution of 



+ 



A 
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Multiplying this equation by (\ -«rUVW*J$)->\ , we obtain 



Thus, 



ui.rA = ( (dom, n (d x ) + A) 1 (fi-in, n (d 1 )f 2 

(d ^Ml.n + /^2,n(^o V2 



Lemma [3.21 ensures that {^A. n {^Q 1 )f2)n weakly converges to rj^{dQ l )f 2 . Thus, by Theorem 
14. 3[ we deduce that (tti, n )n. weakly converges to (do^it^ 1 ) + t4) _1 (/i — )f2) =: Ui- 

Moreover, (dQ 3 Ui >n )n converges pointwise to <5o~ 3 i>i. Using the equality 

U2,n = ~fJ'3,n( S 1 ) U l,n + A»2,n(^0 1 )^o" 1 /a E ^ #2), 

we deduce, with the help of Corollary 13.61 for the first term on the right. hand side and 
Theorem 14.41 for the second term, that 



u 2 ,n - v 2 := -%(^>i + V2(d l )d \f 2 e H u _i(R; H 2 ) 
as n — > oo. We arrive at the limit system 



Vsido 1 ) lj W + VO 0) W V V2(do 1 )do 1 f2 

, we obtain 

^i(^ 1 )+^(5 - 1 )r /2 (a - 1 )- 1 r /3 (a - 1 ) 7/4(0 o - 1 )7/2(0 o - 1 )- 1 \ /V\ , M 0\ A;A 
Next, we consider the operator 



vi(-) + vd-)v2(-y%(-) 7/4(-)%(-) _1 ^ _ A fvi(-) o \ / i o 



r/ 2 (-)- 1 %(-) r/^.)- 1 ; VO 1 A 7/ 2 (-)-V V%(0 1 



By Theorem I4.4[ we deduce that 772O) 1 * s a (c")-material law with strictly positive zeroth 
order term on the range of M 2 2,i(0) for some c" > 0. Moreover, r\\ is a (c')-material law by 
Theorem 14.31 Hence, using Theorem 16.114 we deduce the existence of e", r", d" > such 
that N eH co ' c "'(B(r",r");L(H 1 ®H 2 )) nW°{B{Q,e"); L(H X ®H 2 )). □ 

Remark 4.6. Theorem 14.51 contains a structural result concerning homogenization. We 
have proved that O-analytic material laws lead to O-analytic material laws after the ho- 
mogenization process. Hence, it cannot be expected that the homogenized material law 
contains fractional derivatives with respect to time or explicit delay terms, e.g. a delay 
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4 A general compactness theorem for the homogenization of evolutionary equations 

operator like r_/,/ := /(• — h) for some h > 0, as these are not analytic at 0, see e.g. |2"2"| 
pp. 448 (a),(c)]. By Theorem 12.11 we see that the limit equation is also well-posed and 
causal. The assertion concerning the range of the material law N may be interpreted as 
"the main physical phenomenon remains unchanged under the homogenization process". 
Indeed, the difference between the wave equation and the heat equation written as a first 
order system as in [2"§j Example 1.4.6] or |321 Example 3.2] is the range of the zeroth order 
term in the material law. Due to the limiting process there might be higher order parts in 
the expansion of the material law, which result in memory effects. 

Corollary 4.7. Let H be a separable Hilbert space, e,c, r > 0, A: D(A) ii H — > H skew- 
selfadjoint. Denote by P : H — > N(A) 1 , Q : H — > N(A) the orthogonal projections onto the 
respective spaces N(A) 1 and N(A). Assume that the operator A has the (A r C)-property ; 
i.e., (D(PAP*),\-\pap*) (H,\-\h) is compact. Let (M n ) n be a bounded sequence in 
H m {B(0,e)]L(H))nH m '%B(r 7 r)]L\H)) with M n (0) ^ c on R(M n (0)) = fl(Afi(0)) for all 
nsN. Denote by q 2 : H -» N(Mx(0)) n N(A) 1 , q 4 : H -» 7V(Afi(0)) n N(A) the canonical 
orthogonal projections and assume 

q2M' n (0)q:(q i M^0)qt)- 1 = q^M^)* ^{q A M^)* ^)- x for all n e N. (4.1) 

Then there exists uq > 0, e', d > and (rik)k a strictly monotone sequence of positive 
integers such that for all, v > uq, the sequence 

(d M nk (d^) + A) k 

G- converges to 

d {PMSo^P + P^d^n^rS^P + P*r /4 (5 - 1 )r ?2 (5 - 1 )- 1 g 

+Q*V2(d 1 )-%(do 1 )P + Q*V2{d l y l Q) + A 

in H v -i(R; H), wher% 

!&(■) := lim PM nfc (-)P* - PM nk {-)Q*{QM nk {-Wr l QM nk {-)P* 

>oo 

ife(-) := lim (QM^Q*)- 1 

fe^oo 

%(■) := )wi{QM nh (-)Q*)- 1 QM nti {-)P* and 

k—*co 

774(0 := lim(PM nfc (-)Q*)(gM nfc (.)Q*)- 1 . 

fe^oo 

Proof. The assertion follows by applying Theorem 14.51 to 

PM n (-)P* PM n (-)Q*\\ (PAP*0\\ 
QM n (-)P* QM n (-)Q*J J n >\ OJJ- 

2 The limits are computed in the way similar to Theorem 14 . 5 1 with Hi = N(A) and H2 = N(A). 



{(M n ) n ,A) 



22 



Remark 4.8. The compatibility condition (14.11) may be hard to check in applications. How- 
ever, there are some situations in which the Condition (14. ip is trivially satisfied: 

• A is one-to-one; then N(A) = {0} and = 0. 

• #(Mi(0)) 2 N(A) 1 ; then JV(Mi(0)) n iV(A) 1 = {0} and g 2 = 0. 

• Mi (0) is onto; then the preceding condition is satisfied. We remark here that this con- 
dition was imposed in [29] Theorem 2.3.14]. This condition corresponds to hyperbolic- 
type equations in applications. 

. M;(0) = M;(0)* ; then q 2 M' n {Q)qt{q 4 M' n {Q)qlY l = q 2 M' n ^Y -ql^KW it)' 1 ■ 



We do not yet know whether the compatibility condition is optimal. We can however give 
some examples to show that the other assumptions in (P2) are reasonable. The following 
example shows that without the requirement on A to have the (iVC)-property the limit 
equation can differ from the expressions given in Theorem 14.51 or Corollary 14.71 

Example 4.9 (Compactness assumption does not hold). Let u, e > 0. Consider the 
mapping a : M — > M given by 

a ( x ) := X[o±)( x - k ) + 2 X[i,i](x - k) 

for all x 6 [k, k + 1), where k e 7*. Define the corresponding multiplication operator in 
L2OR), i- e - f° r e C£°(R)j we define a(n ■ m)<j) := (x i-> a(nx)4>(x)) for n 6 N. Note 
that a(x + k) = a(x) for all x e R and k <= Z. Let / e iJ^R; L 2 (R)). We consider the 
evolutionary equation with (M n (5 ( 7 1 )) n := (5 ( T 1 a(n • m)) n and A = i : L 2 (M) — > L 2 (R) : 
(f) h- ► 20. Clearly, iV(A) = {0}. By [3T?] Theorem 4.5] or [TO] Theorem 1.5], we deduce that 



M n ^>(z~z^ a(x) dx^j = (z~^z) eU™(B(0,e);L(L 2 (R))) 



as n — > oo. If the assertion of Corollary 14 . 71 remains true in this case, then (5oM n (5 _1 ) + A) n 
G-converges to | + i. For n 6 N, let u n e H V>0 (M.; L 2 (R)) be the unique solution of the 
equation 

(5 M n (5 - 1 ) + A)u n = {a{n ■ fh) + i)u n = f. (4.2) 
Observe that by [TUl Theorem 1.5] 



= (a(nm) + i) f — - [ \ (a(x) + i) dx ) f =: u. 
as n —*■ oo. We integrate 

(a(ar) + i)" 1 dx = ~(1 + i) -1 + J(2 + 



I 



1 1 



o 2^ ' 2 
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4 A general compactness theorem for the homogenization of evolutionary equations 
Inverting the latter equation yields 

QW) + i)" 1 dx)"' = (1(1 + ,-)- + 1(2 + o-')" 1 - g + 

Hence, it satisfies 

(§+')-'-(£♦&)-'■ 

which of course is a contradiction. 

In the next example, the uniform positive definiteness is violated. 

Example 4.10 (Uniform positive definiteness condition does not hold). Let H = C, v > 

and, for n e N, let M^d^ 1 ) = d Q l \, A = 0, / 6 ^ i0 (M)\{0}. For neN, let it n e P^oQR) 
be defined by 

d M n (do l )u n = -it n = /. 

n 

Then (it n ) n is not relatively weakly compact and contains no weakly convergent subse- 
quence. 

In the final example, the range condition in (P2) is violated. 

Example 4.11 (Range condition does not hold). Let H be an infinite-dimensional, sep- 
arable Hilbert space. Let (4> n ) n be a complete orthonormal system. For n e N define 
M n (5 ~ 1 ) := (0 n , -}</> n + 5 _1 (1 — (0 n , ■)(/)„). For the sequence (M n ) n the range condition 
in (P2) (applied with A = 0) is violated. Let / 6 H ufl (R;H), v > 0. For n e N, let 
it„ g P^qQR; P) be such that 

5 M n (5 ( 7 1 )ii ri = 5 <^ n , u n )4> n + u n - (4> n , u n )4> n = f. 

It is easy to see that {u n ) n is bounded. Take the inner product of the last equation with 
(f) m for some m e N. If n e N is larger than m we arrive at 

(u n , (i) m ) = </, m >, 

and we deduce that (5oM n (5 _1 ))„ G-converges to BqBq 1 = 1. This, however, does not yield 
a differential equation. 
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5 Applications 



We demonstrate the applicability of our main theorem, to the mathematical models of 
some physical phenomena. For notational details, we refer to |29| pp. 34 and p. 98] or to 
[231 3.2 Examples]. 



Thermodynamics 

Let a, 6 M, < a < (3, ft g R N open and bounded, N <=N. Recall [HI Definition 4.11]: 

M(a, /3, ft) := {k e L^fi)***; Re<«;(:r)£, ^ a|£| 2 , ^ /3|£| 2 , f 6 a.e. x 6 ft}. 

Let M se if(a, ft) := {/? 6 M(a, /3, ft); ft selfadjoint a.e.}. For k e M(a,/3,ft) denote 
by «;(m) the associated multiplication operator in L 2 (ft) Ar . Let (« n )„ be a sequence in 
M se if(a;, /3, ft). Recall from Section [T] that a first order formulation of the heat equation 
with Dirichlet boundary conditions in the context of evolutionary equations introduced in 
|21| is the following 

i (A (o o \ / o div\\ f Ul A = (h 

^ I Kntm)- 1 ^ I grad J J U 2 J \f 2 



We want to apply Corollary 14.71 with 
and 



div 
grad 



1 0\ fO 



-i 



The compactness condition on the operator I ^ , ^ V ) has already been established 



,grad 

e.g. in [29| Remark 3.2.2] or |3T| the end of the proof of Theorem 4.3] and the sequence 

(ft n (m) _1 ) n is a sequence of selfadjoint operators. Since M n (0) = Mi(0) = f J for all 

n £ N, the range condition is satisfied. Thus Corollary 14.71 applies. The sequence (M n ) n 
could be replaced by some convolution terms. Moreover it should be noted that the case of 
not necessarily symmetric K n 's has been considered in [31] . There, however, a second-order 
formulation was used which, for more general material laws, may not be available, see also 
the last paragraph in Section [TJ The homogenized equations are derived in Section [JJ see 
equation ( II .Tp . 
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5 Applications 



Electromagnetism 

Let Q E K 3 - The general form of Maxwell's equations in bi-anisotropic dissipative media 
used in [3] is 

dJ\ A + (^^A + f ™*WfE\ (j 



7* n) \cr 2 i* cr 22 *y y— curio J J \H J \0 
with a (c)-material law (see Section [2]for a definition) M(d l ) := ^) °" 12 * 

\3 



7 /iy " \cr 2 i* a 22 * 



for a c > 0. Here, cr^ are L(L 2 valued functions on M, vanishing on IR<o and being such 
that the temporal convolutions <jjk* yield 0-analytic material laws, j, k e {1, 2}. Moreover, 
s 'y \ 

' 1 to col f n *^ t r\inf nn/^ o^T'i 1 1 r nricif n ro rlofimf q in 7^ _^f~)^6 



the operator f # J is selfadjoint and strictly positive definite in L 2 (f2) . We emphasize 

here that the convolution kernels may also take values in the linear operators on L 2 (f2) 3 , 
which are not representable as multiplication operators, thus, in this way, generalizing the 
assumptions in [3]. Now, consider a sequence of (c)-material laws (M n ) n of the above 

form with non-singular, strictly positive zeroth order term: ( £ " ^ n ) ^ d > for all 



ri£N. Then, Corollary 14.71 applies if we assume Q to be bounded and to satisfy suitable 
smoothness assumptions on the boundary, see e.g. |T9 | [23 | 135 ] . Indeed, the range condition 

is satisfied since M n (0) = ( £ ™ ^ n J is onto for all n e N and, since iV(M n (0)) = {0}, 

\7n f^nj 

the compatibility condition fl4.ll) also follows. Note that the homogenized equations are 
more complicated than in the case of the heat equation. This is due to the fact that the 

/ curl\ 

orthogonal projection onto the nullspace of A = I , acts in both components 



curl 

since both the operators curl and curl have an (infinite-dimensional) nullspace. In case of 
the heat equation with Dirichlet boundary conditions the operator grad is one-to-one and 
so the projection onto the nullspace of A in this case only acts in the second component 
and is thus easier to handle. 

To illustrate the versatility and applicability of Theorem 14.51 we show how our methods 
apply to the equations of thermopiezoelectricity. 



Thermopiezoelectricity 

We assume O g R 3 to be open and bounded. The equations of thermopiezoelectricity 
describe the interconnected effects of elasticity, thermodynamics and electro-magnetism. 
The set Q models a body in its non deformed state. We recall the formulation as in |22| 
6.3.3, p. 457], where the model given in |13] is discussed. The unkonwns of the system 
are the time-derivative of the displacement field v, the stress tensor T, the electric and 
magnetic field E and H as well as the temperature distributions 9 and the heat flux Q. 
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Recalling the spatial derivative operators from the introduction and defining Div and Grado 
as the divergence on matrix-valued functions and the symmetrized gradient with Dirichlet 
boundary conditions, respectively, we describe the equations as follows 



ro 



(qq 

C- 1 C- l d C~ X X 
d^C- 1 e + dPC-H p + d*C~ x X 
p 

A*C _1 p* + \*C~ 1 d a + A*C- X A 



V 











\ 



\0 











+ 



q + qi(a + Kd ) 1 J 
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Div 











o\ 


\ 




fv\ 




(f\ 


Grado 





















T 
















curl 












E 




j 
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- curl 















H 






















div 






e 




9 


V o 








o 


grad 




) 








W 



where C, e, n, q , q±, a, k, d, A, p are bounded linear operators in appropriate L 2 (£l)-spaces. 
To frame the latter system into the general context of this exposition, we find that 



.4 



( 


Div 











o\ 


Grado 


























curl 











- 


- curl 


























div 


V o 








o 


grad 


0/ 



The material law is given by 



/eo o 

o c- 1 
o re- 1 
















C~ x d 

e + d*C- x d p + d*C~ x X 
p 
A*C _1 p* + X*C- x d a + A*C _1 A 

yo o ooo 











q + q 1 (a + k8 ) 1 J 

M is a (c)-material law, if we assume that one of the following conditions is satisfied 

(i) £>0i C, e, /x, qo, a — p*e~ 1 p, k are selfadjoint and strictly positive, 

(ii) £oi C, £, p, qi, a — p*s~ 1 p, k are selfadjoint, strictly positive, q\H~ l = K~ l q\ and q = 0. 
Indeed, the material law can be written as a block operator matrix in the form 



m(5 - 1 ; 



Mn 

M 22 (5 - 1 ) 
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6 Auxiliary results on O-analytic material laws 



where 





(Oo 















\ 










o c- 1 


c 


- x d 





c-u 










M n = 


(FC- 1 


£ + d*C~ l 


d 


p + d*C- 


l x 



































\0 \*C~ X 


p* + A*C- 


*d 


a + A*C 


-W 










/l 


0\ 




/ £>0 















/l 


o \ 




1 










c- 1 










1 d 


A 




d* 1 













£ 







OOlOe 


-ip 







1 










fi 







1 







\0 A* 


1 l) 









a- 


p*e 






^0 


1 / 



is strictly positive and, by choosing v > sufficiently large such that 
enough, 



becomes small 



M 2 2(d : ) = q + qi(a + k3 ) 1 = q + qxK 1 d 1 + 



00 

(1=1 



AC a) K 



is such that either M22(0) or M2 2 (0) is strictly positive. Thus, M is a (c)-material law for 
some c > 0. 

Considering a sequence (M n ) n of such material laws with the boundedness and uniform 
positive definitess assumptions from Corollary 14.71 one may derive a homogenization result 
for these equations. We will not do this explicitly here. However, in order to satisfy the 
range condition, one has to assume that all entries of the material law sequence satisfy 
either condition or (jul). To deduce that A has the (iVC)-property, we have to impose 
suitable geometric requirements on Q as in the previous example. 

We refer the intertested reader to more examples of first order formulations of standard 
evolutionary equations in mathematical physics to [211 122] • With these formulations it is 
then straightforward to see when and how our homogenization result applies. 



6 Auxiliary results on O-analytic material laws 

In this section, we provide the remaining theorems needed in Section 0J Our main concern 
will be the discussion of 0-analytic material laws, i.e., material laws that are analytic at 
6 C, cf. [2"T| 3.3 Some special cases]. To establish Theorem 14. 51 similarity transformations 
of 0-analytic material laws have to be discussed, where our main interest is to show that 
under any of these similarity transformations a (c)-material law transforms into a (c')- 
material law for suitable d > 0. In order to achieve the main goal of this section, Theorem 
16.111 some technical results are required. We start with a fact concerning Hardy space 
functions. 
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Lemma 6.1. Let X be a Banach space, e > 0, /i(-) = 2n=o(') ?1 / i n e W°(B(Q, s); X). Then 
for all k,n e Nq we have 

(i) IKII ^ WiAL (f) n 

r«; IK=* ^"Vnii ^ 2 iia»Hoo (?)* allz ^ B (°> t) ■ 

Proof. The first assertion follows immediately from Cauchy's integral formula (integrate 
over a circle around with radius e/2) and the second is a straightforward consequence of 
the first. □ 

With these estimates, we can establish some structural properties of O-analytic material 
laws. Recall that the inner products discussed here are linear in the second and conjugate 
linear in the first component. 

Proposition 6.2. Let H be a Hilbert space, e, c> 0, < r < e/2. Let M be a material law 
m n x (B (0, e)-L(H)) n H x ' c (B(r, r);L(H)). Then M(0) is self adjoint. For cp = fa fa e 
R(M(0))®N(M(0)), the inequalities 

(M(0)fa, fa) ^ and (ReM'(0)fa, fa) ^ fafa, fa) 

/io/<i. //, m addition, R(M(0)) ^ H is closed, there exists d > 0, such that for fa e R(M(0)) 
we have 

(M(0)fa,fa)^d(fa,fa). 

Proof. We expand M into a power series about 0: M(z) = 2^=o z "^n f° r 21 
and suitable (M„) n in L(H). Then M(0) = M and Af'(0) = Mi. For g 
x^ := Im(M(O)0, 0) and := Re(M(O)0, 0). It is easy to see that T : B(r, r) 
iR : z >-* z^ 1 is homeomorphic. Thus, for z £ B(r,r) with Zi := ImT(z), z 2 : = 
we have 

c<0, 0> ^ Re<0, z _1 M (z)^> 

oo 

= Re(zzi + z 2 )(ix 4> + y£ + Re<0, J] z n ~ x M n fa) 

n=\ 

oo 

= + z 2 y<f, + Re<0, ^] z n ~ x M n fa). 

The left-hand side is non-negative. The last term on the right-hand side is bounded for 
z —*■ 0. Moreover, since T is bijective (in particular, for every z 2 the values of z\ range over 
the whole real axis), it follows that x§ = 0. Thus, we arrive at 

X 

c(fa fa < z 2V4> + Re<0, J] z n ~ l M n fa. 

n=l 
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6 Auxiliary results on O-analytic material laws 



Now, since 2 2 can be chosen arbitrarily large, while the second term of the right-hand 
side remains bounded, it follows that y$ ^ 0. Thus, for every (f) e H, we deduce that 
(cp, M(0)fa> ^ 0. Since M(0) is a bounded operator in the complex Hilbert space H, the 
operator M(0) is selfadjoint and positive (semi-)definite and therefore H = R(M(0)) © 
N(M(0)). Let (j) e N(M(0)). Then for e/2 > r] > 0, 

00 

fc=i 

00 

= Re<0,M / (O)0> + r/Re<0,2^" 2 M fe 0>. 

fc=2 

If we let 77 — » 0+, we get that Re<(0, M'(0)<^>) ^ c|0| 2 . Now, M(0) is invariant on its range 
and the restricion of M(0) to its range is one-to-one. Thus, if R(M(0)) is closed, the closed 
graph theorem implies that M(0) : i?(M(0)) — » i?(M(0)) is continuously invertible. By the 
spectral theorem for continuous and selfadjoint operators it follows that M(0) is strictly 
positive on its range. □ 

Remark 6.3. We shall note here that the converse of Proposition 16.21 is also true in the 
following sense: Let M g K x (B(0,e);L(H)) be such that M(0) = M(0)*. Assume that 
there exist d, c> such that for all fa e R(M(0)), fa e N(M(0)) 

(M(0)fa,fa) ^ d(fa,fa) and <ReM'(O)0 2 ,0 2 > ^ c(fa,fa). 

Then R(M(0)) g 5 is closed and for < r ^ — , M e W°' c ^{B(r, r); L(H)), 

2max{^i,d L ] 

cf. [23 Lemma 2.3] or [22J Remark 6.2.7], where ^1 := | (f + f (II^IL f) 2 + f ll M lloc) 
and?:=min{||M||^(|) 2 |,|}. 



Proof. It is easy to see that R(M(0)) is closed. Let (M n ) n in 1/(5) be such that M(z) = 
Iin=o znM n ^ all z e 5(0, e). By Lemma EU we have ||M'(0)|| ^ f and, for all 

< 5 ^ e/A and z e 5(0,5), ||2* =2 ^ n_1 M n || ^ 25 (f ) 2 ||M|| X . For 1/ ^ maxji/!,?- 1 }, 
^ e B(l/(2i/), l/(2z/)), = (fa, fa) g 5(M(0)) N(M(0)) and 77 > 0, 

<0, Rez- 1 M(z)0> 

00 

= (Re^ 1 ) (fa,M(0)fa) + (fa ReM'(0)fa + Re(<j), J] z n ~ l M n fa 

n=2 

^ H0i| 2 + c|0 2 | 2 - 2 ||M'(0)|| \fa\\fa\ ~ ||M'(0)|| |0 X | 2 - 25 \\M\\ X |0| 2 
^ (Vd - 77 ||Af'(0)f - ||M'(0)||) |0 X | 2 + (c - -) |0 2 | 2 - ^|0| 2 
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a ^-,(||M L ?) 2 - M 4-|)w 2 + (|-i)w 2 

If r\ = 3/c, using ^ > 1/1, we obtain 

(Rez -1 M(V)0, 0> ^ T^rf - | (\\M\U - ||ML | _ £ j + £^2 >> £|^». n 

This completes the general discussion on O-analytic material laws. With the main results 
of Section [4] in mind, it is natural to discuss the following situation for material laws M: 

Assumption 6.4. Assume there exist Hilbert spaces Hi,H 2 and constants e, c > 0, < 
r < e/2 with 

M e W°{B{^ e)\ L(Hi H 2 )) n H x ' c (B(r, r);L{H x ® H 2 )) 
and R(M(0)) ^H 1 ®H 2 is closed. 

Before we turn to Gauss-transformations on the material law, we study some structural 
properties of a material law satisfying Assumption 16.41 These properties are stated in the 
next theorem for which we need the following elementary prerequisite. 

Lemma 6.5. Let H\,H 2 be Hilbert spaces. Assume that 

is self adjoint and positive definite. Then M 12 = M 21 and if M 22 = then M 12 = 0. 

Proof. It is easy to see that M n = M* x and M 22 = M 22 and thus 

/ M 12 \ = (Mu M 12 \ _ (M n \ 
\M 21 J \M 21 M 22 ) \ M 22 ) 

is selfadjoint. Assume now that M 22 = 0. If M 12 0, then there exists <f> 2 ) e Hi® H 2 
such that Re(Mi202; <t>\) = (M 12 <j) 2 , 0i) + (M 2 i<j)\, 2 ) < 0. For a > we deduce that 

° S < ffi "o j (at) ' (at) ) " <M "*' *> + " {<M ^ + <AW " ■ 
which yields a contradiction if a is chosen large enough. □ 

In the following, for Hilbert spaces Hi,H 2 we denote the canonical orthogonal projection 
Hi © H 2 —*■ Hj onto the jth coordinate by irj, j e {1,2}. 
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6 Auxiliary results on O-analytic material laws 



Theorem 6.6. Let M satisfy As sumption \6.4\ Using the notation from Assumption \6.4 
we define 

G x := fl(7riM(0)7rf), G 2 := 7V(7nM(0)7r*), G 3 := i?(7T 2 M(0)7r|), G 4 := iV(7r 2 M(0)7r 2 *). 
T/ien M /ias the following form: 



( 



M 



/M[? M x 




(o) ^ 




M 3 ( ? M 3 ( ? 
V 0/ 



+ z 



(Mil\z) M$(z) M$(z) M$(z)\\ 
M&\z) M$(z) Mg\z) M$(z) 
M^{z) M${z) M£\z) M&\z) 

\M£\z) M$(z) M£\z) M&\z)JJ 



6 B(0,e);L 




where for j,k e {1,2,3,4} we have M$ e H cc {B(0 1 e))L(G k ,G j )) and if j,k e {1,3} 



we 



have 



M ( ? e L(G k , Gj). Moreover, there is d > such that Mj 0) ^ d for j e (1, 3}. 



Proof. By Proposition 16. 2[ we know that M(0) is selfadjoint and strictly positive definite 
on its range. Thus, 

is selfadjoint and strictly positive definite and therefore Hi j+ i = Gj Gj+i for j 6 {1, 3}. 

2 "* 2 



We denote by £>j : Hi i7 2 ~* the orthogonal projections onto Gj and define M, 



^M(O)^, M. 



(i) 



(0) 



M 



6 j (M - M(0)) Q* k for all j, fc e (1, 2, 3, 4}. Hence 

+ z 










Ml? 


mJ?\ 










M 2 ( ? 


Mi? 


Mi? 


m£) 


Mi? 
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U? 


Mi? 





o / 



M«(z) M«(*) M«(z)\\ 
M^O) M 2 ( ?(z) m£>(z) m£>(z) 

M 3 W(V) MjJ(z) Mg } (V) M^Oz) 

\m«(«) mI^z) Mi?(z) M«(z)yy 



As M(0) is selfadjoint, Lemma shows that M^?* = Mj? e L(G k ,Gj) for all fc,j e 
{1,2,3,4}. Since M(0) is positive definite, it follows that the block operator matrices 



/mJ? o o a4 0) \ 




Vm|? oo o / 



/0 0\ 

M 2 ( ? 
Mi? Mi? 
\0 0/ 



/0 \ 



M 



(0) 



21 




V0 Mi? 0/ 



(0) 



are positive definite as well. Thus, by Lemma 16. 5[ we deduce that M 14 



M 



(or 



li 



M 



(0) 



23 



M 



(or 



:-!2 



and M 



(o) 



21 



M 



(or 



12 



0. 



o. 
□ 
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For the next theorem we note that, for Hilbert spaces Hi, H 2 and B e L[H 2 , Hi), we have 



1 B 
1 



1 -B 
1 



and 



1 B 
1 



-i 



g a/1 + B + LB 



(6.1) 



Moreover, we need the following lemmas: 



Lemma 6.7. Let H be a Hilbert space. Let T e L(H) be continuously invertible, A,Be 
L{H). If A = T*BT, i.e., A and B are similar, then Re A = T* ReBT and if in addition 
ReB ^ c> then Re A ^ pripr. 

Proof. We have 2ReA = A* + A = T*B*T + T*BT = IT* Re BT . Assume that ReB^c 
for some c > 0. Then, for <p e H, 



(Re Acf), (f) = <Re BT0, T0> ^ c(TcJ), T<f>} ^ 



i2W 



□ 



-we 



Lemma 6.8. Let M satisfy Assumption \6.4\ Using the notation from As sumption \6. 4 
define 

d := R(mM(0)nf), G 2 : = iV(7nM(0)O, G 3 := i?(7r 2 M(0)7r|), G 4 := iV(vr 2 M(0)yr*). 

Lei e L(G 3 ,Gi), nS> e L(G 4 ,Gi), iV 4 ( ? 6 L(G l7 G 4 ), iV|J e L(G 4 ,G 2 ) ; 6 
ft°°(£(0, e); L(G 3 G 4 , Gi G 2 )), A^ e ft°°(B(0, e); L(d G 2) G 3 G 4 )) and 



Then, 



Ni:=\z 



Nv := z 



/V (0) /V (0) \ 



JV- 



13 







+ zN[}\z) 



N (o) N (oy 
iv 41 iv 24 



M:= [ z 



I N fA M{z) ( N ^ z) eH™(B(0,e);L(Hi®H 2 )) 

and R(M(0)) = R(M(Q)), M(0) > d' on its range and ReA^'(O) > c' on the nullspace of 
M(0), where 



d' :=d a 1 + 



N- 



(0) 
13 



+ 



(o) 

13 



and c '.= c I a / 1 + 



A r 



(o) 

24 



+ 



A^ 



(0) 



2i 



with d > fremg i/ie constant of positive definiteness of M(0) on z£s range from Theorem 

E3 
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6 Auxiliary results on O-analytic material laws 



Proof. Using the representation of M given in Theorem 16. 6} we compute Ai(0): 



M(0) 





/ikfff M$ 0\ 


M 3 ( ? Mg^ 

V o o o oj 

(M® Afg? 0\ 


Mi? M<? 

\ o o o o/ 



/ 



V 



N$* 



Nlf 




Hence, A4(0) is similar to a positive definite operator. Moreover, the similarity transfor- 
mation commutes with the projector 



P :-- 



(\ 0\ 

10 

\0 0/ 



Thus, M(0) and A4(0) have the same range and are both strictly positive on it. Indeed, 
M(0) is strictly positive on G\ {0} G% {0} and since the similarity transformation is 
a bijection on d {0} G 3 {0}, M(0) is a bijection on G x {0} G 3 {0} as well. In 
view of Lemma ISTTl and Inequality (16. ip . we deduce A4(0) ^ d! on its range. Next, consider 
(1 — P).A/f(0)(l — P). For this purpose, we compute 
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with suitable N^ k e ^°°(5(0, e);L(G k , G 2 )) (k e {3, 4}) and, similarly, we find 



4,)?) (1 - P) 




1 




\\on. 



(0) = 
24 







1 



+ 



7 




for suitable N^ k2 e ft°°(B(0,e); L(G 2 ,G k )) (k e {3,4}). Note that M'(0) consists only of 
the first order terms of Ai. Since the zeroth order term of M leaves the range of P invariant 



and, due to the structure of the transformations 
most of the terms cancel, we arrive at 



1 

N v {z) 1 



(1-P) and (1-P) 



1 N x {z) 
1 



(1-P)M'(0)(1-P) 
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24 
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N. 



(0) 



24 




1 



(i) 
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(0) M« 
(0) m« 
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(0) M. 
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(0) <(0) M«(0)\ 
(0) M$(0) M«(0) 
M&>(0) 



M. 
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M«(0) Af 43 
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Vo m1 2 1} (0) 



\ 
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,01 
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N, 
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Now, by Lemma [6. 71 and Inequality (16. ip . we see ReA^'(O) ^ d on the nullspace of A4(0). 



□ 



Remark 6.9. Consider the following situation where we apply Lemma 16.81 Let Gj be a 
Hilbert space for j e {1, 2, 3, 4}. Let 

^:-|«-( J f |6?f»(fl(0 >e ) i L(G 1 OG„G,OG 4 )) 
and assume JVg>* = AT 3 ( J } e L(G 3 ,Gi). Moreover, let 

with 4f = iV 3 ( 3 0) e L(G 3 ), e L(G 4 ), iV« E %°°(£(o, e); L(G fc , £?,•)), j, fc e {3,4}. 
Then it is easy to see that N 2 (-)N 3 (-) e H x (B(0, e);L(G 3 ®G A , G 1 ®G 2 )) and N 3 (-)N 2 >(-) e 
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%°°(£?(0, e);L(Gi © G 2 , G 3 © G 4 )) and the following expansions hold 

n,wm = ( n ^ n ^°a + (° i; (o) < 

1 J 1 J V ; l0 7V«(0)iV 4 ( ? 



+ 0(z) 



nWnW q\ ( o 



and 

for — * 0. Now, let M and the Gj's be as in Lemma lBTBl Assume the following compatibility 
condition 

n^nS(o)Y = nS(o)n^. 



Then N\ := N2N3 and Ny := N3N2' satisfy the assumptions from Lemma [6.81 

We now turn to the analysis inverses of material laws. Since we need to estimate the norm 
bounds of these inverses, we observe that, for a Hilbert space H and a continuous linear 
operator B e L(H) satisfying Re B ^ h for some h> 0, B- 1 e L(H). Moreover, using the 
Cauchy-Schwarz inequality, we deduce the estimate 

\\B- l \\^l/h. (6.2) 

Another consequence of Re B ^ h is 

ReS" 1 ^ h/(\\B\\ 2 ). (6-3) 

Further, note that, for appropriate linear operators a,/3,j,5 in Hilbert spaces and A : = 
(a — /35 _1 7) _1 , we have formally 

a /3V 1 ( & -A135- 1 



7 5) ' y-S-^A S^jAfiS- 1 + 8- 1 ) ' 

Lemma 6.10. Let H U H 2 be Hilbert spaces, d,c,e > 0. Let L(H X ) b M$ = M$ * ^ d. 
Moreover, let e W°{B{Q, e); L(H k , Hj)) with ReAf^(O) ^ c. Define 

Then there exists e' > depending on e, WMjW^, , c and d such that, for z e B(0,e')\{0}, 

= ( M 121 (z) -M 121 (z)M^ l (z)Mg\z)- 1 \ 

\-M^{zrHl^{z)M 121 {z) m£\z)-H4{ ) (z)M 121 (z)m£\z)m£\z)-i + z^M% \z)^ ' 
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where 

M 121 (z) := (m£> + z (M${z) - MS(z)M^(z)- 1 M^(z))y 1 
= (m$A 1 + 0{z) and 

On the other hand, for M jk e H x (B(0, e'); L(H k , Hj)) with ReM$(0) ^ c and 



Mj:= \ B(Q,e!)\{Q}3z 
there exists e" > depending on c, d, 



M 



(j, k e {1, 2}) such that, for all z e B(0, e"), 



( MMz) -zM } 2 1 (z)M^{z)Mg\zr 1 \ 

" \-zM£(z)-Hl£{z)M 121 {z) zHdg\z)-H4l\z)M 121 (z)M[l\z)Mg\zr 1 + zMg\z)-\) ' 

where 

M 121 (z) := (Afg> + a (m^z) - Mff^^^)- 1 ^^)))" 1 
= (Aff?)" + 0(z) and 

In particular, M v e H CC {B(Q, e");L(Hi H 2 )). 

Proof. The expressions for the inverses of Mj and Mp can be verified immediately. The 
asymptotic expansions are straightforward applications of the Neumann series expan- 
sion. The respective convergence radii can be estimated in terms of e, ||Mjr|| ,c and d 

or £, Mj k (j, k £ {1, 2}), c and d by Lemma [6TT1 and Inequality (16.21) . □ 

oo 

Theorem 6.11. Let H\,H 2 be seperable Hilbert spaces, c,d,e,r > 0. Let (N n ) n = 
((nH'I En)) be a bounded sequence inH™(B(0,e)\L(Hi®H 2 ))n'H a3 ' c (B(r,r);L(H 1 ® 
H 2 )) and M e V°(£(0, e); L(fTi # 2 )) n H^ c {B{r,r)] L(H l ff 2 )). ^Issnme toa< /or 
all n e N we have R{M(0)) = R(N n (0)) and M(0),N n (0) ^ d on R(M(0)). Denote by 
qj : Hj — > R(ir*) n JV(M(0)) (j 6 {1,2}) £/ie canonical ortho-projections. Assume for all 
fieN £/ie compatibility condition 

(^N^Mq*)- 1 q 2 N^ n (0)qlY = ^ N '^)<& M^JO^)" 1 ■ 
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Then there exist e' , r', d > depending on e, c, r, ((M^ , sup{|| A^H^ ; n e N}, d such that 

(JVi,„)n := ((* -> N 12 , n {z)N 22 , n {z)- 1 )) n 

(N V , n ) n := {{Z ~ iV22 ! n(^)~ 1 A r 21,n(^))) n 

are bounded in / H co (B(0, e')\ L(H 2: Hi)) andW°{B(U,e');L(Hi,H 2 )), respectively. More- 
over, denoting by M 1 e H™{B(0, e'); L{H 2 , #1)) and M v e H™{B(0, e'); L(H U H 2 )) the 
respective limits of (N ltTlk )k and (Nir )nh )^ for a strictly monotone sequence of positive inte- 
gers (rik)k, we have 

1 ±M 1 (z)\ %/rf ^ ( 1 
±M v {z) 1 

e ?T (5(0, e')-L{E x H 2 )) n %™> d {B{r' , r'); L{H X H 2 )), 

and R(M(0)) = R(M(0)). 

Proof. In the following, we use Hilbert spaces Gj, j e {1,2,3,4} as in Theorem 16.61 and 
represent N n for n e N using bounded operators N^ n , (j,k) e {(1, 1), (1, 3), (3, 1), (3, 3)}, 

and Hardy space functions N^ n , j, k e {1, 2, 3, 4}, as in Theorem 16.61 From Lemma 16.101 
we have an explicit expression for N22,n( z )~ 1 , namely 

r (o) N 



w*)- 1 -!^ +ow ,, 0(1) Ifcr^O. 

0(1) ^AfWjO)- 1 + O(l) 

Moreover, we have an estimate for the radius of convergence e' for the Neumann expansion 
involved in this expression in terms of sup{||iV n || ; n e N},d,c,e. In particular, z <— » 

N 22 , n (z) _1 satisfies the assumptions on in Remark 16.91 (note that ^A^^ J = 

{^33 5 by Theorem 16. 6ft . Moreover, using Theorem 16. 6[ we deduce that N\ 2 . n and 
N 2 i } n satisfy the assumptions on N 2 and N 2 > in Remark 16.91 Indeed, we have 

Jr / {N® ()\ ( N$Jz) N$Jz) 

N 12n (z) := 0w JVl 3,« +2 )?\ n ; |f\ nV ; 

e ^(^(O, e); L(G 3 G 4 , Gi © G 2 )) 
,'iVj? 0\ ( N$Jz) N$Jz 



e 7*°°(fl(0, e); L(Gi G 2 , G 3 G 4 )) 



with ^A^ n ^ = ATg^ by Theorem 16.61 Moreover, Remark 16.91 shows that Ni >n and Ny tn 
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satisfy the assumptions imposed on Ni and Ny in Lemma ESI More precisely, we have the 
expansions 

o Jv£>(o)jy£>(o)-' + () 



r'«* "».« " | +o(.-) 



for suitable continuous linear operators Ni^ n , A^ 41n . We deduce that 

Moreover, the compatibility condition is precisely 

<.(0)<„(0)- lN )* = < n (0)- 1 iVl^ n (0). 



Lemma 13.31 together with the fact that computing the adjoint is a continuous process in 
the weak operator topology ensures that Mi and My satisfy the assumptions imposed 
on N\ and Ny in Lemma 16.81 To estimate the norm bounds of M\ and My in terms 
of d, c, sup{ || A^H^ ; n e N} and e, we use Lemma [6.11 and Lemma 16.101 Hence with the 
help of Lemma I6.8[ we may estimate the constants of positive definiteness of M. (0) and 
Re A^'(0) on R(M(0)) and 7V(.M(0)), respectively, also in terms of d, c, sup{|| A^IL ; n e N} 
and e. Note that we also have i?(.M(0)) = R(M(0)). Remark 16.31 implies the remaining 
assertion. □ 
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